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Abstract 

The  feasibility  of  a  re-entry  guidance  scheme  that  utilizes  ballistic 
flight  to  attain  a  variable  surface  range  is  investigated.  The  equations 
of  motion  are  derived,  and  the  re-entry  trajectory  is  divided  into  four 
phases:  pull  up,  skipout,  ballistic,  and  terminal.  The  pull  up,  skipouf, 
and  terainal  phases  are  inside  the  atmosphere,  while  the  ballistic  phase 
is  outside  the  atmosphere.  The  conjugate  gradient  optimization  method  is 
used  in  the  skipout  phase  to  determine  the  re-entry  vehicle  control  needed 
to  achieve  the  desired  initial  ballistic  phase  conditions.  The  cost  function 
uses  the  errors  between  the  actual  trajectory  and  a  reference  trajectory 
which  is  obtained  by  projecting  a  ballistic  path  into  the  atmosphere.  The 
equations  for  the  reference  trajectory  are  presented.  Ballistic  theory  is 
discussed  and  a  parametric  selection  technique  i3  developed  as  a  means  of 
determining  the  ballistic  parameters  needed  for  a  specific  range.  An  op¬ 
timal  control  solution  is  found  for  one  set  of  ballistic  initial  conditions. 
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2.  Introduction 


Background 

Consider  a  spacecraft  orbiting  the  noon,  prior  to  a  return  flight  to 
Earth.  The  time  of  departure  from  lunar  orbit  is  dependent  upon  many 
variables,  one  of  which  is  the  Earth's  rotation.  That  is,  the  location  of 
the  landing  site  moves  with  the  Earth.  The  spacecraft  must  enter  the  at¬ 
mosphere  within  a  certain  time  period  so  that  the  landing  site  will  be 
within  range  of  its  guidance  capabilities.  In  the  case  of  the  recent 
Apollo  missions,  the  maximum  attainable  surface  range  was  approximately 
22 GT;  nautical  miles  (nr.)  (Ref  1).  Alternate  landing  sites  were  available 
for  the  Apollo  splashdowns,  but  every  alternate  site  requires  the  use  of 
additional  recovery  ships  and  other  resources. 

To  eliminate  the  need  for  excessive  backup  landing  areas,  it  is  de¬ 
sirable  to  have  a  guidance  scheme  that  will  allow  a  re-entry  vehicle  (RV) 
more  flexibility  in  surface  range.  Ideally,  a  spacecraft  returning  from 
the  soon  would  be  able  to  land  at  one  fixed  location,  independent  of  the 
distance  involved,  including  circumnavigation  of  the  glebe. 

One  method  of  attaining  more  surface  range  is  by  including  a  segment 
of  ballistic  flight  xn  the  re-entry  trajectory.  If  the  spacecraft  is  made 
to  skipout  c*  the  atmosphere  after  initial  entry,  it  would  then  be  on  a 
ballistic  trajectory.  Assuming  the  craft  re-enters  the  atmosphere  due  to 
the  Earth’s  gravitational  influence,  the  surface  distance  traveled  during 
ballistic  flight  would  depend  only  upon  the  state  of  the  vehicle  at  atmos¬ 
pheric  exit  (skipout  conditions).  The  skipout  conditions,  in  turn,  would 
be  a  function  of  the  total  range  desired.  flfnen  the  term  range  is  used 
alone,  it  will  naan  surface  range). 
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Purpose  and  Scope 

The  purpose  of  this  thesis  Is  to  investigate  the  feasibility  of  a 
re-entry  guidance  scheme  that  will  utilize  ballistic  flight  to  achieve  a 
variable  surface  range.  The  surface  range  refers  to  the  total  distance 
between  the  initial  point  of  re-entry  and  the  target  point  above  the 
landing  site. 

The  Bain  emphasis  of  this  study  is  placed  upon  1)  finding  the  sklpouc 
conditions  for  a  desired  ballistic  trajectory,  and  2)  optimally  driving 
the  states  of  the  RV  tc  those  skipout  conditions.  The  control  necessary 
to  obtain  the  correct  states  is  found  by  an  optimisation  scheme. 

The  only  type  of  re-entry  considered  directly  in  this  study  is  one 
resulting  from  a  lunar  return  flight,  but  extensions  to  other  types  could 
be  made.  Also,  the  g  forces  and  the  vehicle  nose  heating  are  not  directly 
constrained. 

Method  of  Analysis 

The  problem  was  attacked  by  dividing  the  trajectory  into  four  dis¬ 
tinct  phases.  The  theory  of  ballistic  flight  was  then  considered  in  order 
to  establish  the  necessary  skipout  conditions.  It  was  decided  that  the 
normal  method  of  calculating  ballistic  range  'see  equation  (10))  would 
not  suffice  for  the  problem  at  hand,  because  fcu?  range  was  dependent  upon 
two  arbitrary  variables.  Consequently,  a  linear  stepping  technique  was 
devised  which  would  allow  ballistic  range,  and  hence  the  skipout  conditions, 
to  be  an  indirect  function  of  only  &  stepping  variable.  Then,  to  simplify 
the  analysis,  the  terminal  phase  was  arbitrarily  assigned  s  fixed  value 
cf  range.  This  was  done  by  assuming  that  a  terminal  phase  guidance  scheme 
could  land  the  RV  for  a  realistic  range  of  state  conditions  occurring 
after  ballistic  flight. 
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Emphasis  was  then  placed  upon  finding  a  net  hod  (guidance  scheme)  to 
attain  the  various  skipout  conditions.  Several  unsuccessful  methods  were 
initially  tried.  The  first  method  was  an  attempt  to  derive  a  control  law 
from  the  derivative  of  Q.  Equation  (7)  was  solved  for  the  angle  of  attack 
and  dQ/ds  was  defined  e a  the  slope  of  a  straight  line  between  the  actual 
Q  and  the  desired  Q  at  skipout .  This  control  law  failed  because  it  could 
not  control  the  initial  p.'-ft  of  the  flight,  due  to  the  amount  of  energy 
involved. 

The  second  idea  was  to  hold  the  derivative  of  the  flight  path  angle 
at  zero  when  the  RV,  during  its  upward  flight,  acquired  the  needed  skipout 
angle.  The  idea  was  only  partly  successful  because  once  the  angle  began  to 
stray  froa  the  desired  value,  the  errors  could  not  be  corrected.  In  an 
effort  to  make  this  second  idea  self-corrective,  it  was  converted  into 
a  second  order  approximation  by  taking  the  second  derivative  of  &  with 
respect  to  range.  The  resulting  coefficient  of  d5/ds  was  set  equal  to 
2?^  and  a  control  law  was  found  by  solving  for  the  angle  of  attack  (which 
was  contained  in  the  coefficient).  The  control  law  operated  as  a  bang- 
bang  controller  in  an  attempt  to  function  properly.  Since  a  violent  con¬ 
trol  is  undesirable,  this  method  was  put  aside. 

The  conjugate  gradient  optimization  technique  was  finally  chosen  as 
a  suitable  means  of  finding  a  control  that  would  drive  the  RV  to  the  proper 
conditions.  The  conjugate  gradient  cost  function  for  this  problem  Is 
based  upon  the  differences  between  the  actual  states  of  a  re-entry  vehicle 
and  the  states  of  an  imaginary  ballistic  vehicle  whi« h  is  traveling  on  a 
ballistic  path  extended  into  the  atmosphere.  For  a  guidance  scheme  to 
function  properly,  the  actual  trajectory  should  converge  to  the  projected 
path  before  ballistic  flight  begins. 


Overview 

The  mathematical  issodels  are  set  up  in  Chapter  II.  The  items  presented 
are  the  simplifying  assumptions,  the  atmospheric  model,  the  vehicle  lift 
and  drag  equations,  a  summary  of  the  equations  of  motion  (or  states)  and 
the  vehicle  heating  and  deceleration  equations. 

Chapter  III  then  goes  into  the  problem  formulation.  A  typical  re¬ 
entry  trajectory  is  divided  into  phases  and  a  method  of  range  allocation 
for  the  ballistic  phase  is  developed.  Included  in  the  chapter  is  dis¬ 
cussion  of  gone  ballistic  equations,  development  of  a  reference  trajectory, 
derivation  of  the  equations  needed  ia  the  conjugate  gradient  method,  and 
discussion  of  the  computer  requirements. 

Chapter  IV  then  covers  the  results  of  this  study,  while  Chapter  V 
contains  the  conclusions. 

There  are  three  appendices  to  this  report.  The  equations  of  motion 
are  derived  in  Appendix  A  and  the  conjugate  gradient  theory  is  given  in 
Appendix  B,  Appendix  C  shows  a  plot  of  angle  of  attack  versus  the  lift 
coefficient. 
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II.  Mathematical  Model 


The  assumptions  listed  below  were  made  in  order  to  reduce  the  com¬ 
plexity  of  the  problem  while  retaining  the  most  important  characteristics. 

1.  All  motion  is  two-dimensional. 

2.  The  Earth  is  non-rotating  and  is  spherical  in  shape. 

3.  The  only  bodies  considered  are  the  Earth  and  the  re-entry  vehicle. 

4.  The  Earth's  atmosphere  above  an  altitude  of  400,000  feet  is 
neglected. 

5.  Below  400,000  feet  in  altitude,  the  atmospheric  density  is  con¬ 
sidered  a  function  of  altitude  only,  as  expressed  by  the  model 


P  -  Pft  e 


where  pQ  is  sea  level  density 


pc  <*  0.0026703  slug3/ft^ 


and  8  is  an  atmospheric  constant 


B  -  O.OOC 0425211877  ft' 


Aerodynamic  Model 

In  addition,  it  is  assumed  that  control  of  the  re-entry  vehicle 
will  be  accomplished  by  directly  changing  the  angle  of  attack  of  the 
vehicle.  The  resulting  lift  and  drag  specific  forces  can  be  represented 


L  -  I  i  p  V2  Ct 
2  a  L 


■IS 


M 
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D  -  I  i  p  V2  Cn 
2  a  v 


O) 


S  ? 

The  inverse  wing  loading,  “,  Is  1.18070  ft  /slug, 
efficients,  for  a  hypersonic  lifting  body  (Ref  2), 
Hewtoaiaa— Flat-Plate  drag  polars  and  are  as  shown. 


the  lift  and  drag  co~ 
are  represented  as 


*  1.82  sina  cosa  jsinaj 
Cjj  -  0.042  +  1.46  |sin3a| 


(4) 

(5) 


For  this  vehicle  model  then,  (L/D) _  -  2 

’  max 

The  assumptions  and  models  are  used  in  the  derivation  of  the  equa¬ 
tions  of  motion. 


Equations  of  Motion 

The  equations  of  motion  were  written  relative  to  an  inertial.  Earth 
centered,  coordinate  system.  The  equations  (see  Appendix  A  for  the  de¬ 
rivation)  are  represented  as  derivatives  with  respect  to  surface  range,  s, 
and  they  are  listed  here  for  continuity. 


dxj  ^  dh 
ds  ds 


—  tan  5 


(6) 


»  (tan6  (2— Q)  -  Rp  £  seed  0  CD)  ~  (7) 

ds  ds  m  Re 


~3-  -  ii6_  -  (  I  _  l  -  rp  S _  gcc6  Clj  1_  (8) 

ds  ds  Q  2a  Re 
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where  h  is  the  altitude  above  sea  level.  Q  is  an  energy  paraaeter  (defined 
in  Chapter  III),  6  is  the  flight  depression  angle,  and  R  **  Re  +  h. 

Re  «  20925738  ft 

is  the  radius  cf  the  Earth,  taken  from  Bate  (Ref  3). 

Heating  and  Deceleration  Equations 

Although  vehicle  nose  heating  and  deceleration  forces  were  not  directly 
constrained,  they  ware  observed  to  insure  that  the  magnitudes  were  reason-' 
able. 

The  total  heat  was  obtained  by  integration  of  tne  heat  rate 

—  "  C.  Pi  V2  seed  —  (8a) 

ds  ueat  Rg 

7  1 

2  7  7 

where  Cbeat  »  2  x  10~®  BTU*sec  /ft  ‘slugs 
is  the  heating  coefficient. 

Deceleration  was  computed  frora  the  time  derivative  of  velocity,  i.e., 

sin5  -  D  (62) 

dt  R2 

where  \i  »  1.407654  x  10^  ft^/see2  (froa  Ref  3)  is  the  Earth's  gravita¬ 
tional  constant. 
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III.  Problem  Formulation 


Phase  Definition 

To  simplify  the  problem  analysis,  the  re-er.try  trajectory  is  separ¬ 
ated  into  four  distinct  phases  which  are  shown  in  Fig.  1  and  are  defined 
below.  A  flat  Earth  representation  is  used  in  the  figure  for  ease  of 
illustration. 

1.  Pull  up  Phase  -  -  That  part  of  the  trajectory  from  initial  re¬ 
entry  at  400,000  feet  to  the  point  where  the  flight  path  angle,  6,  ini¬ 
tially  changes  sign  (from  positive  to  negative). 

2.  Skipout  Phase  -  -  That  part  of  the  trajectory  from  the  point 
where  5  changes  sign,  as  mentioned  above,  to  the  point  where  the  RV  leaves 
the  atmosphere.  The  guidance  scheme  is  used  in  this  phase. 

3.  Ballistic  Phasa  -  That  part  of  the  trajectory  occurring  above 

400,000  feet,  wnere  ballistic  theory  governs  flight.  The  state  values  at 
the  beginning  of  this  phase  are  referred  to  as  the  skipout  (so),  or  skip, 
conditions. 

4.  Terminal  Phase  -  -  That  part  of  the  trajectory  following  ballis¬ 
tic  flight.  It  begins  at  400,000  feet  and  ends  with  a  landing  on  Earth. 

Ballistic  Theory 

Ballistic,  or  free  flight  theory  is  not  new.  Equations  describing 
ballistic  motion  indicate  that  ranges  up  to  halfway  around  the  globe 
(approximately  10,800  no)  are  possible  for  values  of  Qso  <_  1.  For  Qgo  >  1, 
greater  ranges  are  possible  (see  Ref  3). 

Q  is  a  non-flioensional  parameter  defined  as 


trajectory  Phase  Definition 
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The  values  of  Qg0  and  680  (skipout  conditions)  are  used  to  predict  the 
range  attainable  for  a  particular  ballistic  flight.  Since  ©any  different 
valued  of  range  are  possible  for  the  many  different  combinations  of  Qao 
and  6S0,  the  ballistic  phase  is  where  most  nf  the  range  adjustments  will 
be  made. 

An  important  feature  of  a  ballistic  trajectory  is  that  it  is  symmetric 
(see  Fig.  5).  If  Q  and  5  are  known  at  the  beginning  of  free  flight*  the 
reflected  skipout  conditions  will  exist  at  the  end  of  free  flight,  pro¬ 
vided  the  initial  and  final  radii  are  equal.  The  vehicle  in  this  study 
will  leave  and  enter  the  atmosphere  at  the  same  altitude,  so  the  states 
of  the  RV  at  the  end  of  the  skipout  phase  must  be  a  reflection  of  the 
states  at  the  beginning  of  the  terminal  phase. 


Range  Allocation 

It  was  mentioned  in  the  previous  section  that  variation  of  the  skip¬ 
out  conditions  causes  a  range  variation  between  the  beginning  and  end 
points  of  ballistic  flight.  The  free  flight  range  angle,  $■,  is  normally 
calculated  from 


arccos 


1  -  Q: 


SO 


COS 


2  6 


so 


/l  +  Q__(QC_  -  2)cos2  6 

so  ^so  so 


(10) 


However,  equation  (10)  depends  upon  two  variables,  Qso  and  6co,  both  of 
which  may  be  arbitrarily  varied.  Also,  there  is  no  unique  combination  of 
Qs0  and  fiso  for  many  choices  of  desired  range  angle. 

Finding  a  unique  solution  is  possible  by,  in  effect,  reducing  the 
number  of  variables  to  one.  This  is  accomplished  by  using  a  linear 
stepping  technique,  which  was  arbitrarily  selected  in  the  form  shown  be¬ 
low. 


r. 


o 
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QSo  *  Ql  +  (Q2  -  Ql>A 


6eo  m  61  +  (62  *  61)& 


or  aa  alternate  to  equation  (12)  could  be  written 


«9  -  («,  -  6.)A 


-so  2  "  *°2  "  °1 


where  A  is  a  parametric  stepping  variable  that  ranges  from  0  to  1,  and 
the  numerical  subscripts  refer  to  the  upper  and  lover  bounds  cm  Qso  and 
^so* 

After  selecting  upper  and  lower  bounds,  the  values  of  and  6  „ 

’  xGO  BO 

necessary  for  a  desired  range  angle  are  dependent  only  upon  the  value  of 


Equations  (11),  (12),  and  (13)  are  plotted  in  Pig.  2  for  an  arbitrary 


choice  of  values,  i.e. 


■  -9  £  <!„„  £  -  <*5 


S1  -  2‘  £  «SO  £  r 


Notiee  that  only  the  magnitude  of  the  flight  depression  angle  is  needed 
for  this  discussion.  In  Pig.  2,  tne  symbol  'x*  corresponds  to  equation 
(12),  the' symbol  ’O'  corresponds  to  equation  (13),  and  the  small  symbol 
’A1  corresponds  to  equation  (11). 

Using  the  stepping  variable  to  find  values  of  Qso  and  6go,  from 
equations  (11)  and  (12)  or  (13),  corresponding  values  of  range  angle  are 


IK 


%  V. 

•fl 

‘•I; 

If 


2^ 


\ 
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-Pj 

•J 
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generated  froa  equation  (10).  The  range  angles,  as  a  function  of  the 
stepping  variable,  are  plotted  in  Fig.  3.  In  that  figure,  the  symbol  *x’ 
is  associated  with  the  'x'  in  Fig.  2,  and  si&dliarly  for  the  other  axis. 
Referring  to  Fig.  3  and  choosing  either  *x’  (equation  (12))  or  ’0*  (equa¬ 
tion  (13))  for  a  particular  range  depends  upon  the  range  desired  and  the 
slope  of  the  range  curves. 

For  short  ranges,  equation  (12)  would  be  used  to  find  6so,  and  for 
long  ranges,  equation  (13)  would  be  used.  For  any  intermediate  range,  the 
range  curve  with  the  least  slope  would  dictate  which  equation  would  provide 
the  least  sensitivity  to  terminal  errors  in  dgo.  Therefore,  a  logical 
division  of  ranges  is  shown  in  Fig.  4,  which  is  merely  F'g.  3  with  certain 
segments  deleted.  Given  a  desired  rang;.,  Fig.  4  indicates  which  equation 
is  needed  for  calculation  of  the  flight  path  angle  since  the  symbol  *x* 
implies  use  of  equation  (il)  and  the  symbol  ’0*  implies  use  of  equation 
(12). 

By  examining  Fig.  4,  and  the  computer  printout  of  data  (not  shown), 
one  can  find  that  the  ballistic  range  for  the  current  choices  of  Qj_,  Q2, 

&l ,  and  62  is  34,  S®  (2070  nm)  to  346*  (20,760  nm). 

The  total  range  angle,  o,  is  dependent  upon  how  much  distance  the 
remaining  three  phases  can  add.  For  simplicity,  it  would  be  desirable 
for  one  or  more  of  the  ether  phases  to  have  a  fixed  value  of  attainable 
range. 

The  terminal  phase  is  singled  out  for  this  purpose,  because  the  only 
item  of  interest  here  is  how  much  range  can  the  terminal  guidance  scheme 
cover.  For  the  purpose  of  this  study,  it  is  assumed  that  a  guidance 
scheme  exists  which  will  provide  a  nominal  range  of  10C0  no  (chosen 
arbitrarily).  The  assumed  guidance  must  be  able  to  use,  an  initial  conditions. 


Ballistic  Phase  Bangs  Division 
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the  values  of  Q  and  f  vhich  exist  at  the  beginning  of  the  terminal  phase 
due  to  free  flight  symmetry. 

The  range  contributed  by  the  pull  up  phase  and  the  skipout  phase  are 
determined  later  from  computer  data. 

Reference  Trajectory 

The  cost  function  chosen  for  the  optimization  method  is  based  on 
errors  between  the  actual  trajectory  and  a  reference  trajectory.  (The 
coat  function  is  presented  in  the  next  section). 

The  purpose  of  this  section  is  to  present  the  equations  used  to  ob¬ 
tain  a  reference  trajectory  for  the  skipout  phase  of  flight.  This  refer¬ 
ence  path  is  sisply  defined  as  the  projection,  into  the  atmosphere,  of  the 
desired  ballistic  phase  trajectory  (see  Fig.  5). 

To  project  the  ballistic  flight,  one  can  use  classical,  two-body, 
orbital  mechanics  (Ref  3)  and  the  fact  that  the  ballistic  trajectory  forms 
part  of  an  orbit.  First,  the  values  that  remain  constant  for  a  particular 
orbit  (orbital  constants)  must  be  found.  The  orbital  constants  are: 
specific  mechanical  energy,  E*, 

E*  -  -  £  (14) 

2  R 

specific  angular  momentum,  h*, 

h*  *  R  V  cos5  (15) 

the  eccentricity,  e*, 

- 

V  y<* 

and  the  length  of  the  serai-latus  rectum,  p*, 

p*  «  **£ 
y 


(16) 


(17) 
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A  Ballistic  Trajectory  (Froa  3) 
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The  angle  fro®  periapsis  to  the  position  radius  vector  is  v. 


f 

* 

v  “  arccos 

i-  £l-l 

e*  l  a 

• 

<• 

If  the  values  of  Vso,  Rso,  6so  are  used  in  equations  (14),  (15),  and  (18), 
the  values  of  £*,  h*,  e*,  p*,  and  vso  are  then  known  and  can  be  used  as 
shown  next. 

Let  be  the  total  range  angle  covered  by  the  skipout  phase.  Let 
a 2  be  the  actual  range  angle  traveled  by  the  RV  at  an  '  point  in  the  skip- 
out  phase  (see  Fig.  6).  Then,  as  the  RV  travels  along,  it  has  a  particu¬ 
lar  value  for  R,  Q,  and  6  for  each  value  of  02 .  For  the  same  value  of 
02,  a  corresponding  R,  Q,  and  6  can  be  found  for  an  imaginary  vehicle 
traveling  along  Che  projected  ballistic  trajectory.  The  equations  used 
tw  calculate  the  states  of  the  imaginary  vehicle  are  given  below. 


vb  *  vso  -  < °1  -  °2> 


(19) 


Rb  - 


1  +  e*  cos  V}, 


(20) 


vb  -  -\p > 


(21) 


% 


?b*b 


(22) 


5v  »  arccos 


h* 


*b  vb 


(23) 
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Actual  and  Reference  Trajectory  During  Skipout 
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The  differences  between  the  states  of  the  two  vehicles  (one  imaginary) 
are  the  errors  used  in  the  cost  function. 


Conjugate  Gradient  Equations 

The  equations  necessary  to  apply  tha  conjugate  gradient  method  to 
the  skipout  phase  are  derived  in  this  section,  using  the  theory  which  is 
explained  in  Appendix  B. 

The  cost  function  to  be  minimized  is 


J 


8f  _ 

.5  /  [e?Ae  +  (Act)B(Aa)  ]ds  + 


(.5  ?p  el 


S*Se 


(24) 


The  term  added  to  the  integral  is  i  penalty  function  on  the  terminal 
conditions.  The  symbol  An  is  rhe  change  in  control,  and  the  vector  e  con¬ 
tains  the  errors  between  the  states  of  the  actual  and  the  reference  tra¬ 
jectory. 


- 

Xj^s)  ■  hb^s) 

e  2 

eQ 

xs 

x2(s)  -  Q{j(b) 

e6 

«* 

*3(5)  -  Ajj(s) 

(25) 


The  weighting  matrices  are  chosen  to  be 


Ax  0  0 

0  A2  0 

0  0  A, 


(2b) 


B  «  scalar 


(27) 
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P  - 


0 

0 


0  0 
f2  0 

0  P3 

^  4* 


(28) 


Thus,  the  cost  can  be  written  as 


minimize  J  *>  .5  /  (e^  A2  +  A^  +  (Aa)2  B)ds 


(29) 


+  .5<e*  Fj  +  F2  +  e*  ?3>6.S£ 


To  convert  to  a  Meyer  formulation,  the  integrand  cf  the  integral  portion 
of  J  is  made  into  a  state  and  its  differential  equation  is 

dx*  2  2  2  2 

-  ■  *5'eh  Aj  +  eg  A2  +  ej  A3  +  (Ao)  B)  (30) 

ds 


The  other  states  are  repeated  here  for  reference. 


d*l 

ds 


R 


tan$ 


Re 


(6) 


dx2 

ds 


(tanfi  (2-0)  -  Ro 


—  seed  Q 
a 


“D- 


1 


(7) 


^*1  »  (  I  -  1  -  Rp  —  sec6  C»  )  — —  (8) 

ds  Q  2a  Re 


Forming  the  haoiltonian  from  equations  (75),  (30),  (6),  (7),  and  (8) 
yields 


H  "  -  Ai  —  tan 5 
1  *e 
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+  X4  .5(e2  Ax  +  e2  A2  +  e2  A3  +  (Aa)2  B) 

The  adjoints,  from  equations  (76)  and  (31),  are 

dli  S 

— —  “  [X^  tan6  +  P  g  Q  sec6  (^(I-RB) 

dn 

+  X3  p  1-  seed  Cl(1-RS)]  I-  -  X4  eh  Aj 
2a  R« 


(31) 


(32) 


* 


! 

! 

t 

I 

J 


o 


3“-  -  [X2(tanfi  +  Rp  |  seed  C^)  +  X3  --  ]  ~ 


-  X4  eQ  A2 


(33) 


dXi  o  o  S 

gg-  •  [X^  R  sec^d  -  X2(sec  6(2-Q)  +  Rp  ^  Q  Cjj  tand  seed) 


o 


+  X3  *>  li  CL  tan6  sec{5  I"  -  x4  eS  A3 


dX^ 

ds 


(34) 

(35) 


X^  *  constant 

The  transversality  conditions  are  obtained  by  using  equation  (77), 
where 

$(x(sf))  -  x4(sf)  +  .5(e2  T1  +  *2  F2  +  e2  *3}^  (36) 
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therefore 


h(sf>  ’  <s 

(37) 

W  ■  (cQ  Vs-a, 

(38) 

W  *  <*5 

(39) 

*4(sf)  -  1 

(40) 

Equations  (35)  and  (40)  imply  that  «  1,  therefore,  no  integration 

of  equation  (35)  is  required. 

Finally,  equations  (78)  and  (31)  are  used  to  find  an  expression  for 
the  gradient,  which  oust  equal  zero  for  an  optimal  solution. 


G(r)  •  -  (8-76  \2  Q  sino  cosa  +  1.82  A3(2cos2a-sin2a)) 
*  (seco  sina  sgn(a)  p  |_  Rj  (Aa)B 


(41) 


The  equations  derived  in  this  section  were  used  in  the  conjugate  gradient 

algorithm  (Appendix  B)  to  obtain  an  optimal  solution  for  the  skipout  phase 
of  flight. 


Computer  Methods 

The  necessary  computer  programs  were  written  in  Fortran  IV  language 
and  were  executed  on  a  CDC  6600  computer. 

The  pull  up  and  skipout  phases  were  programmed  separately.  For  the 
pull  up  phase,  the  states  were  integrated  forward  until  the  flight  path 
angle  went  to  zero  degrees.  The  states  at  that  point  were  used  as  the 
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initial  conditions  for  the  skipout  phase, 

For  the  skipout  phase,  the  conjugate  gradient  algorithm  shown  in 
Appendix  B  was  used.  The  state  and  adjoint  equations  were  integrated  using 
a  fixed  step,  Runge-Kutta  method.  A  fixed  number  of  points  were  necessary 
to  score  the  state  values,  which  were  subsequently  used  during  reverse 
integration  of  the  adjoint  equations.  Good  results  were  obtained  using 
arrays  of  801  points  each  for  the  states,  ad joints,  gradient,  direction 
of  search,  control,  and  perturbated  control  (used  in  the  alpha  search). 

The  total  conjugate  gradient  program  used  42K  of  core  memory. 

For  a  set  of  skip  conditions,  the  penalty  terms  were  initially  set 
to  zero  while  the  integral  term  of  the  cost  function  (see  equation  (29)) 
was  minimized.  Then  the  penalty  terms  were  applied. 

The  conjugate  gradient  algorithm  was  terminated  when  the  cost  did 
not  vary  in  six  significant  digits. 
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IV.  Results 

The  re-entry  conditions  of  Apollo  10  ware  used  as  the  initial  condi¬ 
tions  for  the  pull  up  phase.  A  constant  angle  of  attack  (control)  of 
54.74®  was  applied  to  drive  the  flight  depression  angle  to  zero.  This  value 
of  control  produces  the  maximum  coefficient  of  lift  (see  Appendix  C). 
Consequently,  the  RV  does  not  penetrate  the  atmosphere  too  deeply  and  the 
nose  heat  rate  peak  tends  to  be  minimized.  The  pull  up  phase  end  condi¬ 
tions  are  summarized  in  Table  I,  where  the  states  are  h,  0,  and  5.  The 
value  for  velocity  was  calculated  from  Q,  and  the  total  nose  heat  is  pro¬ 
vided  as  a  check  for  those  interested  in  physical  limitations.  The  largest 
deceleration  force  in  the  pull  up  phase  was  approximately  4  g’s,  which 
occurred  as  the  flight  path  angle  vent  to  zero. 

The  pull  up  phase  final  conditions  shown  in  Table  I  were  used  as 
initial  conditions  for  the  skipout,  phase.  The  skipout  phase  final  condi¬ 
tions  were  found  by  using  a  desired  ballistic  range  a-.d  the  parametric  range 
selection  equations,  as  explained  in  Chapter  III.  To  make  the  values  of 
Qso  and  5g0  come  out  nice  the  desired  ballistic  range  was  chosen  as  } 6,924 
nm  or  282.07  deg,  therefore,  from  Figures  4  and  2,  ’RANGE  ANGLE  l’  and 
’FLT  PATH  ANGLE  1*  were  used.  For  the  above  range,  the  stepping  variable  is 
0.4  which  gives  skipout  conditions  or  Qgo  **  1.1,  6go  »  -4  deg,  and  of 
course  hgo  *  400,000  ft. 

The  conjugate  gradient  technique  was  then  applied  for  the  selected 
end  conditions,  using  an  initial  guess  of  control  of  71*  (constant).  The 
usual  difficulties  were  encountered  in  application  of  the  optimal  theory 
and  in  finding  suitable  values  for  the  weighting  matrices  in  the  cost 
function,  equation  (29).  The  problem  with  the  matrices  was  that  they 
were  sensitive  to  the  skipout  conditions  and  the  skipout  range  selected. 
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Table  I 

Pull  up  Phase  End  Conditions 

(values  rounded) 

Variable 

Initial 

Condition 

Final 

Condition 

Units 

B 

0.00 

380. 0C 

nst 

h 

400,000.00 

218,259.00 

ft 

Q 

7.00 

1.72 

- 

(5 

6.61 

-0.00 

deg 

V 

36,303.00 

33,871.00 

ft/sec 

<1 

0.00 

11,178.00 

BTU/ft2 

Table  II 

Cost  Function  Weighting  Values 


Variable 


Value 


.1  x  10“10 
,1  x  101 
.1  x  103 
.1  x  10° 

.1  x  103 
,1  x  1013 
.1  x  1015 
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Here  specifically,  a  sec  of  weightings  that  minimized  the  cost  under  one 
set  of  terminal  conditions,  would  not  work  using  another  set,  i.e.,  the 
alpha  search  could  find  no  minimum.  The  weightings  finally  selected  are 
shown  in  Table  II,  and  after  a  few  trial  runs,  a  skipout  range  of  600  am 
was  finally  used. 

Convergence  to  the  skipout  conditions  required  40  iterations,  which 
took  600  seconds  of  execution  time.  The  conditions  were  reached  with  less 
than  1%  error,  which  resulted  in  a  ballistic  phase  range  error  of  r.pprox- 
imately  +80  na.  The  convergence  errors  are  summarized  in  Table  III,  and 
the  optimal  states  are  shown  converging  on  the  reference  trajectories  in 
Figures  7,  8,  and  9.  The  optimal  control  is  presented  in  Fig.  10. 

The  maximum  deceleration  in  the  skipout  phase  was  8  g's  and  the  peak 
heat  rate  was  385  BTU/ft^*sec.  Both  of  these  occurred  when  the  control 
went  to  90*,  which  gives  maximum  drag  (see  Fig.  10),  at  the  beginning  of 
the  trajectory.  The  total  nose  heat  at  the  end  of  the  optimal  skipout 
trajectory  was  approximately  24,400  ETU/ft  .  The  skipout  phase  time  of 
flight  was  132  seconds. 

The  approximate  surface  range  traveled  for  the  entire  flight  is 
totaled  in  Table  IV. 
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Optimal  Altitude  Trajectory 


Optimal  Q  Trajectory 


Optimal  Flight  Path  Angle  Trajectory 


fig.  10.  Optimal  Control  Trajectory 


Range  (ran) 

Range  (deg) 

Pull  up  Phase 

380 

6.4 

Sklpout  Phase 

600 

10,0 

Ballistic  Phase 

17,000 

283.3 

Terminal  Phase 

1,000 

16.6 

Totals 

18,980 

316.3 
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V.  Conclusions 

A  t/s-entry  guidance  scheme  that  utilises  ballistic  flight  to  attain 
more  surface  range  is  feasible,  although  the  deceleration  forces  and  the 
peak  heat  ra*e  may  be  high  at  the  beginning  of  the  skipout  phase  due  to 
the  sudden  application  of  full  drag. 

The  total  surface  range  available  depends  upon  the  values  of  Qso  and 
6go  selected  for  the  ballistic  phase  and  upon  the  surface  range  assigned 
to  the  other  three  phases.  The  minimum  and  maximum  ballistic  range  is 
calculated  in  Chapter  Ill  for  0.9  £  Qso  £  1.4  and  2*  <_  6go  _<  7*.  Using 
the  other  ranges  listed  in  Table  IV,  the  total  surface  range  available 
for  re-entry  is  between  67.5*  (4050  nm)  and  379*  (22,740  nm).  The  max¬ 
imum  range  is  greater  than  one  circumference  of  the  Earth.  For  different 
boundaries  on  Qgo  and  6so  in  the  range  selection  equations,  shorter  or 
longer  total  distances  are  possible. 

The  conjugate  gradient  technique  is  suitable  for  finding  an  optimal 
control,  but  the  weighting  values  in  the  cost  function,  equation  (29), 
change  as  a  function  of  the  skipout  phase  terminal  conditions.  Consequent¬ 
ly,  a  significant  change  in  terminal  conditions  requires  changes  in  the 
weighting  value*. 

However,  due  to  the  long  execution  time  involved  for  convergence  (600 
seconds),  the  conjugate  gradient  method  is  an  unrealistic  choice  as  an 
onboard,  real  time  controller.  A  real  tine  controller  is  needed  which  is 
capable  of  finding  the  control  necessary  for  the  skipout  phase.  The  compu¬ 
tation  should  be  performed  during  pull  up  and  could  be  based  on  expected 
pull  up  phase  terminal  conditions.  Adjustments  would  have  to  be  made  for 
any  errors  between  the  expected  terminal  conditions  and  the  actual  terminal 


conditions. 
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Any  further  investigations  in  this  area  should  be  aade  using  an 
ticial  scheme  that  is  better  suited  than  the  conjugate  gradient  method 
meet  terminal  conditions. 
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Appendix  A 

Derivation  of  State  Equations 

The  equations  that  describe  the  motion  of  an  unpowered  lifting 
body  (re-entry  vehicle)  in  the  Earth's  atmosphere  are  derived  in  this 
appendix,  using  the  assumptions  made  in  chapter  II. 


Coordinate  System 

The  geometric  relations  used  are  shown  in  Fig.  11.  The  coordinate 
frame  E  is  a  locally  inertial.  Earth  centered  frame  where  axis  is 
aligned  with  the  initial  point  of  atmospheric  re-entry. 

The  coordinate  frames,  in  addition  to  being  right-handed,  orthogonal, 
and  cartesian,  were  chosen  such  that  each  vector  would  be  aligne.  with 
a  certain  axis.  The  coordinate  frames  rotate  with  their  respective 
vectors,  therefore,  they  also  serve  as  a  reference  for  the  angle  of 
attack  a(t)  and  the  flight  depression  angle  6(t). 

From  Fig.  11  the  vectors  can  be  expressed  in  matrix  form  as 


(radius) 

(velocity) 

(lift) 

(drag) 

(gravity) 


(42) 

(43) 

(44) 

(45) 

(46) 
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e2 


Earth 
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where  the  subscript  refers  to  the  coordinate  frame  in  which  the  vector 
is  measured. 

The  angular  rotation  rates  are 


or,  in  skew-symmetric  form 


M 


a 


0  -  Ac 

dt 

£*£  0 
dt 


0 


(47) 


(48) 


(49) 


N 

’m 


U). 


o  - 


dS 

dt 


dS 

dt 


0 

0 


0 

0 

0 


(50) 


where  the  equation  is  the  angular  velocity  of  the  M  frame  with  respect 
to  the  E  frame. 

The  coordinate  transformations  are 


<5 


,N 

'M 


cos  a 
sine 
0 

cosd 

sin<5 

0 


-sino 
cos  a 
0 

■sini 
cos  6 
0 


0 

d 

0 

0 

0 

1 


(51) 


(52) 
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where  transformation  is  from  the  frame  denoted  by  the  superscript  to  the 
frame  denoted  by  the  subscript. 


Derivation  of  the  Time  Derivatives 

The  next  step  is  to  equate  velocity  (acceleration)  components  ob¬ 
tained  through  coordinate  transformations*  with  velocity  (acceleration) 
components  obtained  by  using  the  theorem  of  Coriolis  (Ref  4:  Chap  2). 
Therefore,  the  desired  equations  are 


VM  ’  VH 

(53) 

VH  '  +  “H  *!( 

(54) 

*H  ■  °ii  %  +  Hi  +  “h 

(55) 

*K  *  +  “K  VH 

(56) 

Using  the  appropriate  forms  of  equations  (42),  (43),  (47),  and  (52)  to 
substitute  into  the  right-hand  side  of  equations  (53)  and  (54),  snd 
equating  the  results,  one  obtains 


-V  sin  6 

r  •  ~\ 
R 

V  cos  6 

m 

Ro 

0 

0 

- 

From  the  first  row  of  equation  (57) 

R  «  —  =  -7  sin  6 
dt 


Likewise,  the  second  row  yields 


a 


do 

dt 


—  cos  6 
R 


(57) 


(58) 


(59) 


Following  the  same  procedure  for  equations  (55)  and  (56),  one  obtains 
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"o‘ 

-  “2  cos  5  +  L  +  V  (o  +  6)" 

V 

m 

^7  sin  6  —  D 

R^ 

0 

0 

Ose  of  equation  (59)  and  rearrangement  of  row  1  provides 

t  di  u  .  V  t 

6  *»  —  *  -  cos  5  -  _  cos  i  -  ~ 

it  VS2  R  V 


The  final  equation,  which  is  given  in  row  2  of  equation  (60),  is 


(60) 


(61) 


dV 

dt 


~  sin  6  -  D 
R2 


(62) 


The  equations  just  derived  represent  the  motion  of  a  re-entry 
vehicle  with  respect  to  time.  The  equations  will  subsequently  be  re¬ 
ferred  to  as  the  state  equations  or  states,  and  they  are  summarized 
below. 


-  -V  sin  6 
dt 


(58) 


—  ®  X  cos  6 
dt  R 


—  “  iL_  sin  5  -  D 
dt  R2 


(59) 

(62) 


— -  ■  -JL-  cos  6  -  X  cos  6  -  It  (61) 

dt  VR2  R  V 

where,  in  equation  (58),  it  was  decided  to  use  altitude  rather  than 

radius.  Since  R  *  Rg  +  h,  the  time  rate  of  change  is  the  same  for  either 
variable. 


Change  of  Independent  Variable 

Time  is  the  independent  variable  in  the  states  above,  but  in  the 
problem  formulation,  surface  range,  s,  is  of  interest.  Conversion  from 
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time  to  range  as  the  independent  variable  is  accomplished  «sing  equation 
(59)  and  the  fact  that 

a  «  o  Se  (63) 

Also,  the  range  angle,  a,  in  equation  (59)  does  not  appear  in  any  other 
state  equation,  so  it  can  be  eliminated  as  a  state. 

The  inverse  of  equation  (59)  is 

dt  .  _ — E -  (64) 

do  V  cos  5 

or 

dt  - - E - do  (65) 

V  cos  6 

The  derivative  of  equation  (63)  is,  after  rearrangement 

do  ■  ~  (66) 

Re 

Substituting  for  do,  equation  (65)  becomes 


dt  - - S - ds  (67) 

V  cos  6  Re 

Replacing  dt  in  equations  (53),  (61),  and  (62),  and  after  multiplying, 
the  new  states  become 

db.  »  -  E_  tan  5  {  6  ) 

ds  Re 

dV  -  (  _JL.  tan  6  -  -O -  )  I_  (68) 

ds  R  V  V  cos  6  Ra 

d£  ,  (  _JL.  -  1  -  J LI _  )  i_  (69) 

ds  V2R  V2  cos  6  Re 

with  surface  range  as  the  independent  variable.  Since  the  time  of  flight 
ia  not  of  primary  interest  here,  the  equations  are  now  in  a  more  useful 
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fora.  Also,  the  number  of  equations  to  be  integrated  has  been  reduced 
from  four  to  three.  However,  one  change  remains  to  be  cade. 

Q  as  a  State 

The  states  ere  currently  written  in  terms  of  R,  V,  and  6,  but  the 
charts  and  equations  which  cover  ballistic  theory  are  written  in  terms 
of  Q  and  6  (see  Chap  III).  Consequently,  it  is  advantageous  to  have  Q 
as  a  state  rather  than  V.  The  derivation  of  dO/ds  follows,  starting  with 
the  definition  of  Q. 

Q  2  V2R/*i  (9  ) 

The  derivative  with  respect  to  range  is 

«  I  (V2  ®  +  2S  ?  )  (70) 

ds  it  ds  ds 

Substituting  equations  (6  )  and  (68)  for  the  derivatives  on  the  right 
side,  a ad  rearranging,  yields  the  final  form. 


iP  *  (  ten  i  (2 ~Q)  ~  Sp  S.  sec  6  Q  CD)  •*—  (7  ) 

as  ®  Re 


the  value  of  velocity,  if  needed,  is  found  from  equation  (9),  i.e. 


V  *  /Cu/R 


(71) 


Equations  <f>8)  and  (6?)  must  now  be-  rewritten  in  terais  <*e  t’.v-  new 
state,  Q.  The  final  fora  of  the  slat  k  equations  used  for  simulation  are 
eisaEftsrised  here  e-ad  in  Chapter  II. 


4x,  db  8 

— •-£.  *  , —  «r  «  —  tar.  o 

ds  ds 

*  (tan  5  (2-Q)  -  Splnvi  t)  C^I  i- 
do  da  m  R 

e 


(6  ) 

) 
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Appendix  3 

Conjugate  Gradient  Theory 

The  general  theory  of  the  conjugate  gradient  minimization  technique 
is  covered  in  this  appendix.  Also,  the  computational  algorithm  used  t:o 
find  the  optimal  control  is  presented.  Hucls  of  the  material  below  is 
-'‘ran  directly  from  Ref  5. 

Gtaffd 

the  co-'j-fg-it-s  gradient  "wthc-i  5}  is  a  mini aizat ion  technique 

that  can  he  applied  to  optimal  control  problems.  A  '•.lotion  requires 
knowledge  of  the  gradient  trajectory,  its  norm,  and  the  actual  direction 
of  search.  The  search  direction  is  found  by  use  of  the  norm  to  modify 
the  gradient  direction,  consequently,  each  search  is  in  the  conjugate 
direction.  Thus,  even  for  a  poor  initial  guess  of  the  optimal  control, 
the  method  tends  to  converge. 

There  are  a  couple  of  disadvantages  to  the  conjugate  gradient  tech¬ 
nique.  First,  the  method  applies  only  to  unconstrained  problems.  How¬ 
ever,  if  terminal  conditions  or  inequality  constraints  are  present,  pen¬ 
alty  functions  can  be  used  to  convert  the  problem  to  an  unconstrained 
fora.  A  second  disadvantage  is  that  it  is  unable  to  distinguish  between 
a  local  and  a  global  minimum,  which  is  a  fault  common  to  most  minimiza¬ 
tion  schemes. 

The  general  problem  formulation  is  for  the  Meyer  fora  of  an  optimal 
control  problem. 


minimize  J  « 

$(x(tf>) 

(72) 

subject  to  x  * 

f (x,u,t) 

(73) 

x(t0)  » 

constant 

(74) 
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where  x  Is  an  n-order  state  vector,  u  ia  an  s-order  control  vector,  and 
the  initial  time  tQ  and  the  final  time  tjt  are  fixed.  Subsequent  discussion 
of  u  will  consider  only  a*) . 
the  Hamiltonian  is 


H-  *  hh 

i-1 


<?5) 


where  the  ad joints  are 


-  -  t  A«  ,  i-1,. 


d\< 

- -  -  -  i  A« 

at  j»:  axi 


(76) 


and 


Al(tf)  -  ~~  |  *  i»l*« 

3x£  t*tf 


(77) 


The  gradient  is 

G(u)  -  ||  (78) 

A  necessary  condition  for  x  to  be  an  optimal  state  trajectory  is 
that  the  gradient  equal  zero  for  the  optimal  control  trajectory.  Thus, 
the  optimal  control  minimizes  the  Hamiltonian.  The  algorithm  described 
in  the  next  section  is  an  iterative  process  used  to  find  that  control. 


Computational  Algorithm 

After  the  Hamiltonian,  adjoint,  and  gradient  equations  have  been 
derived,  the  optimal  control  caa  be  found  by  following  the  steps  listed 
below.  The  subscript  i  refers  to  the  previous  value. 

1.  Select  an  arbitrary  control  trajectory  (often  chosen  as  a  con¬ 
stant)  and  use  it  to  integrate  the  state  equations  forward.  Proceed  to 
step  (3). 

2.  Perform  an  alpha  search.  This  consists  of  finding  a  distance 
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o' ,  such  chat 

J<ui+1)  "  J(ui  +  a*  8i>  (79> 

is  ainiei2ed.  This  requires  forward  integration  of  the  states  until 
equation  (79)  is  satisfied. 

3.  Integrate  the  adjoint  equations  backwards  as  a  function  of  the 
final  state  values. 

4.  Compute  the  gradient  from  the  states  and  adjoints. 

5.  Calculate  a  new  direction  of  search  s’  according  to 


s*  «  ~G  +  B’  s’ 
i+1  <+l  i+1  i 


(80) 


with 


where  the  norm  is  defined  as 


(81) 


!  |G|H  -  ftfG*  dt  (32) 

co 

For  the  _.iitlal  iteration,  s^  ® 

6.  Repeat  steps  (2-5)  until  the  value  of  cost  in  step  (2)  changes 
less  than  soma  tolerance.  The  current  value  of  control  ir  step  (2)  will 
be  the  optiaal  control. 
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